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We characterize and classify quantum correlations in two- 
fermion systems having 2K single-particle states. For pure 
states we introduce the Slater decomposition and rank (in 
analogy to Schmidt decomposition and rank), i.e. we de- 
compose the state into a combination of elementary Slater 
determinants formed by pairs of mutually orthogonal single- 
particle states. Mixed states can be characterized by their 
Slater number which is the minimal Slater rank required to 
generate them. For K — 2 we give a necessary and suffi- 
cient conditions for a state to have a Slater number one. We 
introduce a correlation measure for mixed states which can 
be evaluated analytically for K — 2. For higher K, we pro- 
vide a method of constructing and optimizing Slater number 
witnesses, i.e. operators that detect Slater numbers for some 
states. 
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I. INTRODUCTION 



In recent years a lot of effort in Quantum Infor- 
mation Theory (QIT) has been devoted to the charac- 
terization of entanglement, which is one of the key fea- 
tures of quantum mechanics ||. The resources needed 
to implement a particular protocol of quantum informa- 
tion processing (see e.g. [0) are closely linked to the en- 
tanglement properties of the states used in the protocol. 
In particular^ entanglement lies at the heart of quantum 
computing |S|]. The most fundamental question with re- 
gard to entanglement is: given a state of a multiparty 
system, is it entangled or not (i.e. is it separable ||)? 
If the answer is positive, then the next question is how 
strong the entanglement is. For pure states in bipartite 
systems the latter question can be answered by looking 
at the Schmidt decomposition ||, the i.e. decomposition 
of the vector in a product basis of the Hilbcrt space with 
a minimal number of terms. For mixed states already the 
first question is notoriously hard to answer. There exist, 
however, many separability criteria, such as the Peres- 
Horodecki criterion [0,9, and more recent concepts such 
as entanglement witnesses and the corresponding "entan- 
glement revealing" positive maps p|JlO||- 

While entanglement plays an essential role in quan- 
tum communication between parties separated by macro- 
scopic distances, the characterization of quantum corre- 



lations at short distances is also an open problem, which 
has received much less attention so far. In this case 
the indistinguishable character of the particles involved 
(electrons, photons,...) has to be taken into account. In 
his classic book, Peres || discussed the entanglement in 
elementary states of indistinguishable particles. These 
are symmetrized or antisymmetrized product states for 
bosons and fermions, respectively. It is easy to see that 
all such states of two-fermion systems, and as well all 
such states formed by two non-collinear single-particle 
states in two-boson systems, are necessarily entangled in 
the usual sense. However, in the case of particles far 
apart from each other, this type of entanglement is not 
of physical relevance: "No quantum prediction, referring 
to an atom located in our laboratory, is affected by the 
mere presence of similar atoms in remote parts of the 
universe" ||. This kind of entanglement between indis- 
tinguishable particles being far apart from each other is 
not the subject of this paper. Our aim here is rather 
to classify and characterize the quantum correlations be- 
tween indistinguishable particles (in our case fermions) 
at short distances. We discuss below why this problem 
is relevant for quantum information processing in various 
physical systems. 

For indistinguishable particles a pure quantum state 
must be formulated in terms of Slater determinants or 
Slater permanents for fermions and bosons, respectively. 
Generically, a Slater determinant contains correlations 
due to the exchange statistics of the indistinguishable 
fermions. As the simplest possible example consider a 
wavefunction of two (spinless) fermions, 
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with two orthonormalized single-particle wavefunctions 
4>(r), x( 7? )- Operator matrix elements between such sin- 
gle Slater determinants contain terms due to the antisym- 
metrization of coordinates ("exchange contributions" in 
the language of Hartree-Fock theory). However, if the 
moduli of cj)(r), xif) have only vanishingly small over- 
lap, these exchange correlations will also tend to zero 
for any physically meaningful operator. This situation is 
generically realized if the supports of the single-particle 
wavefunctions are essentially centered around locations 
being sufficiently apart from each other, or the particles 
are separated by a sufficiently large energy barrier. In 
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this case the antisymmetrization present in Eq. (|l|) has 
no physical effect. 

Such observations clearly justify the treatment of indis- 
tinguishable particles separated by macroscopic distances 
as effectively distinguishable objects. So far, research in 
Quantum Information Theory has concentrated on this 
case, where the exchange statistics of particles forming 
quantum registers could be neglected, or was not speci- 
fied at all. 

The situation is different if the particles constituting, 
say, qubits are close together and possibly coupled in 
some computational process. This the case for all pro- 
posals of quantum information processing based on quan- 
tum dots technology Jll]-[l3|] . Here qubits are realized by 
the spins of electrons living in a system of quantum dots. 
The electrons have the possibility of tunneling eventu- 
ally from one dot to the other with a probability which 
can be modified by varying external parameters such as 
gate voltages and magnetic field. In such a situation the 
fcrmionic statistics of electrons is clearly essential. 

Additional correlations in many-fermion-systems arise 
if more than one Slater determinant is involved, i.e. if 
there is no single-particle basis such that a given state 
of N indistinguishable fermions can be represented as 
an elementary Slater determinant (i.e. fully antisymmet- 
ric combination of N orthogonal single-particle states). 
These correlations are the analog of quantum entangle- 
ment in separated systems and are essential for quantum 
information processing in non-separated systems. 

As an example consider a "swap" process exchanging 
the spin states of electrons on coupled quantum dot s by 
gating the tunneling amplitude between them pgJDjj . 
Before the gate is turned on, the two electrons in the 
neighboring quantum dots are in a state represented by 
a simple Slater determinant, and can be regarded as dis- 
tinguishable since they are separated by a large energy 
barrier. When the barrier is lowered, more complex cor- 
relations between the electrons due to the dynamics arise. 
Interestingly, as shown in Refs. [ fH^Ulljfl , during such a 
process the system must necessarily enter a highly corre- 
lated state that cannot be represented by a single Slater 
determinant. The final state of the gate operation, how- 
ever, is, similarly as the initial one, essentially given by 
a single Slater determinant. Moreover, by adjusting the 
gating time appropriately one can also perform a "square 
root of a swap" which turns a single Slater determinant 
into a "maximally" correlated state in much the same 
way In the end of such a process the electrons can 
again be viewed as effectively distinguishable, but are in 
a maximally entangled state in the usual sense of distin- 
guishable separated particles. In this sense the highly 
correlated intermediate state can be viewed as a resource 
for the production of entangled states. 

We expect that similar scenarios apply to other 
schemes of quantum information processing that involve 
cold particles (bosons or fermions) interacting at micro- 
scopic distances at which the quantum statistics becomes 
essential. For instance, it should be of relevance for quan- 



tum computing models employing ultracold atoms in op- 
tical lattice s |14| , or ultracold atoms in arrays of optical 
microtraps Jl 5[| . 

It is the purpose of the present paper to analyse the 
above type of quantum correlations between indistin- 
guishable fermions in more detail. However, to avoid 
confusion with the existing literature we shall reserve in 
the following the term "entanglement" for separated sys- 
tems and characterize the analogous quantum correlation 
phenomenon in non-separated fermionic systems by the 
notions of Slater rank and Slater number to be defined 
below. 

We are going to formulate analogies with the theory 
of ent anglem ent, and translate several very recent re- 
sults |lO| , |l6|Jl7| concerning standard systems of distin- 
guishable parties (Alice^ Bob) to the case of indistin- 
guishable fermions. In general we will deal with a system 
of two fermions each of which live in a 2if-dimensional 
single-particle space. 

The plan of the paper is as follows: In Section [n] we 
discuss pure states, and formulate the analog of Schmidt 
decomposition and rank - Slater decomposition and rank. 
We then discuss a simple operational criterion for the 
case of two electrons in two neighboring quantum dots 
(if = 2) to determine whether a given state is of Slater 
rank 1. This criterion was first derived in Ref. Il3(|. In 



Section III we define the concept of Slater number for 
mixed states. For K = 2 we present necessary and sufi- 
cient condition for a mixed state to have the Slater num- 
ber 1. This is an analog of the Peres-Horodecki crite- 
rion |?],U in the Wootters formulation . In Section IV 
we extend the results of Section III and define a Slater 
correlation measure which is the analog of the entangle- 
ment formation measure ]l{J. This quantity can be cal- 
culated analytically for the case K — 2, in analogy to the 
Wootters result Jig] , In Section V we turn to the case 
K > 2 and introduce Slater number witnesses of canon- 
ical form (defined in analogy to entanglement |]]r| and 
Schmidt number [^(||l7j witnesses). We construct exam- 
ples of such /c-Slater witnesses, which provide necessary 
conditions for a given state to have the Slater number 
smaller than fc; we also discuss optimization of Slater 
witnesses. Finally, we analyze the associated |2l]] positive 
maps. We close by discussing further analogies, but also 
differences, between entanglement in separated systems 
of distinguishable particles as opposed to quantum cor- 
relations in non-separated systems of indistinguishable 
particles. 



II. SLATER RANK OF PURE STATES 

We consider two indistinguishable fermions each of 
which lives in the single-particle Hubert space C 2K . This 
situation is given, e. g., in a system of two electrons in K 
neighboring quantum dots where only the orbital ground 
state of each dot is taken into account. Alternatively one 
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may think of, say, two quantum dots with an appropriate 
number of orbital states available for the two fermions. 

The states (density matrices) in such a system are pos- 
itive self-adjoint operators acting on the antisymmetric 
space A(C 2K <g> C 2K ). Let us first consider pure states, 
i.e. projectors on a vector e A(C 2K <g> C 2K ). Let 
fai fai a = 1) • • • ) 2if j denote the fermionic annihilation 
and creation operators of single-particle states forming an 
orthonormal basis in C 2K , and |fi) denotes the vacuum 
state. Each vector in the two-electron space can be rep- 
resented as |*) = J2 a ,b w abfift\&), where w ab = -w ba is 
an antisymmetric matrix. We have the following gener- 
alization of the Theorem 4.3.15 from Ref. p^ |, which will 
allow us to define the fermionic analog of the Schmidt 
decomposition: 

Lemma 1 For any antisymmetric N x N matrix A ^ 
there exists a unitary transformation U' such that A — 
U' ZU' T , where the matrix Z has blocks on the diagonal, 



Z = diag [Z ,Zi,. . . , Z M ] , Z = 0, Zi 







Zi 





(2) 



and Z is a (TV - 2M) x (N - 2M) null matrix. 



Proof: Let A be a TV x AT, complex, antisymmetric matrix 
acting on C , A = —A T , hence A? = — A* . Let us define: 
B := AA* = -AAL B is hermitian, B = B\ hence 
diagonalizable by a unitary transformation: B — UDW , 
UW = 1, D - diagonal. Now consider C := WAU*. It 
is easy to check that C is antisymmetric, C T = —C, and 
normal CC^ = C^C. Let us decompose C into its real 
and imaginary parts: C = F + iG; F,G are real TV x N 
matrices. Since C is antisymmetric, so are F and G. 
Since C is normal, F and G commute. Thus F and G are 
real, antisymmetric, commuting matrices. Hence they 
can be simultaneously brought to block diagonal forms 
by a real orthogonal transformation |2^| , F = OFmO 7 ' , 
G = OG bd O T , O is a TV x N matrix, OO t = I, where 



F bd = diag [X , Xi, . . . , X K ] 
G 6d = diag [y ,*i, ■••,*£], 



(3) 



and JVo, Yq are null matrices of some dimensions, Xq = 0, 
Yq = 0, whereas Xj, 1^ are standard antisymmetric 2x2 
blocks: 



Xi 







Vi 


—Xi 




-Vi 






(4) 



Thus C = OZO T where Z has the form (§) and, finally 
A = UCU T = UOZO T U T = (UO)Z(UO) T = U'ZU' T 
with U' unitary. □ 

Lemma 2 Every vector in the antisymmetric space 
A(C 2K ® C 2K ) can be represented in an appropriately 
chosen basis in C 2K in a form of the Slater decomposi- 
tion 
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, TV, form 



where the states f , 

an orthonormal basis in C 2K , i.e. each of these single- 
particle states occurs only in one term in the summation 
(||). The number of nonvanishing coefficients 2, (i.e. the 
number of elementary Slater determinants required to 
construct \^)) is called the Slater rank. 

Proof: Let |f) = J2a,b w abflf b \ty ■ Note that the change 
of basis in C 2K corresponds to a unitary transformation 
of fermionic operators, /J = ^ b U ba (f' b )\ which implies 
that in the new basis w' = UwU T . From Lemma 1 we 
may choose U such that w' will have the form (|^), which 
provides the Slater decomposition. □ 

From the point of view of applications in quantum dot 
computers, it is important to be able to distinguish states 
with Slater rank 1 (which can be easily prepared and de- 
tected) from those that involve more than one elementary 
Slater determinant. In general, given |^) in some basis, 
in order to check the Slater rank, one has to per form 
the Slater decompositon. As we know from Ref. |L3|, the 
situation is simpler for the case K = 2, where we have: 



Lemma 3 [Ref. (T3J] A vector |* 
in A(C 4 O C 4 ) has the Slater rank 1 iff 



Ea.l 



=1 w ab flfl\n) 



T,(|*)) = 



t abCd WabW cd 

a,b,c,d 



= o, 



(6) 



where ( abcd denotes the totally antisymmetric tensor in 
C 4 ®C 4 ®C 4 ®C 4 . 

Remark: The quantity r,(|*)) can be constructed from 
the dual state 



I*) = X> o6 /t/+|n) 



defined by the dual matrix 



IV 



cd ' 



With these definitions we have 



7,(|*)) 



(7) 



(8) 



(9) 



The proof of this Lemma was presented first in Ref. [^3| . 
An alternative proof can be given using Lemma 1 and 
observing that 



det w 



(10) 
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where w is the antisymmetric 4 x 4-matrix defining | \& ) . 

In the Appendix we list some further useful proper- 
ties of ^(l^ 1 )) and the relation of the dualization oper- 
ation to an antiunitary implementation of particlc-hole- 
transformation. An interestring further question are pos- 
sible generalizations of the above result to the case of K 
fermions having a single particle space C 2K . 

III. SLATER NUMBER OF MIXED STATES 

Let us now generalize the concepts introduced above 
to the case of mixed states. To this end, we define the 
Slater number of a mixed state, in analogy to the Schmidt 
number for the case of distinguishable parties : 

Definition 1 Consider a density matrix p of a two 
fermion system, and all its possible convex decomposi- 
tions in terms of pure states, i.e. p = Y^iPil^T)^^]} 
where denotes the Slater rank of the Slater num- 

ber of p, fc, is defined as fc = min{r max }, where r max is 
the maximum Slater rank within a decomposition, and 
the minimum is taken over all decompositions. 

In other words, k is the minimal Slater rank of the pure 
states that are needed in order to construct p, and there 
is a construction of p that uses pure states with Slater 
rank not exceeding fc. 

Many of the results concerning Schmidt numbers can 
be transferred directly to the Slater number. For in- 
stance, let us denote the whole space of density matrices 
in A(C 2K <g) C 2K ) by Six, and the set of density matri- 
ces that have Slater number fc or less, by Sl k . Sl k is 
a convex compact subset of SIk) a state from Sl k will 
be called a state of (Slater) class fc. Sets of increas- 
ing Slater number are embedded into each other, i.e. 
Sl\ C Sl-2 C ...Slk--- C SIk- In particular, Sl\ is the 
set of states that can be written as a convex combination 
of elementary Slater determinants; SI2 is the set of states 
of Slater number 2, i.e. those that require at least one 
pure state of Slater rank 2 for their formation, etc. 

The determination of the Slater number of a given state 
is in general a very difficult task. Similarly, however, as 
in the case of separability of mixed states of two qubits 
(i.e. states in C 2 ® C 2 ), and one qubit-one qutrit (i.e. 
states in C 2 Cg>C 3 ) pj, the situation is particularly simple 
in the case of small K. For K = 1 there exists only one 
state (a singlet). For K = 2 we will present below a nec- 
essary and sufficient condition for a given mixed state to 
have a Slater number of one. One should note, however, 
that in the considered case of fermionic states there ex- 
ists no simple analogy of the partial transposition, which 
is essential for the theory of entangled states. In fact, the 
Peres-Horodecki criterion in 2 x 2 and 2x3 spaces 
says that a state is separable iff its partial transpose is 
positive. It is known, however, that the Peres-Horodecki 
criterion is equivalent to Wootters' result |L8|], relating 
separability to a quantity called concurrence, which is 



related to eigenvalues of a certain matrix. This latter 
approach can be used to characterize fermionic states in 
A(C 4 ® C 4 ). We have the following Theorem: 

Theorem 1 Let the mixed state acting in A(C A (8 C 4 ) 
have a spectral decomposition p — Yll=i l^ r i)(^ r i|) where 
r is the rank of p, and the eigenvectors \^>i) belonging to 
nonzero eigenvalues Xi are normalized as (Wjl^j) = XiSij. 

Let l^i) = J2 a b w abflfb\ty m some basis, and define the 
complex symmetric r x r matrix C by 

C ij = J2e abcd wi b wi d , (11) 

abed 

which can be represented using a unitary matrix as 
C = UCdU T , with Cd = diag[ci, C2, . . . , c r ] diagonal and 
I cx I > |ca| > ... > \c r \. The state p has Slater number 1 
iff 

ki|<f>|. ( 12 ) 

i=2 

Proof: Let us assume that a state p acting in A(C A ®C 4 ) 
has Slater number 1, i.e. 

r r 

p = £|**K*<I = X>*><M (is) 

i=l fc=l 

where all have Slater rank 1, whereas r' can be an 
arbitrary integer > r. But \4>k) can be represented as 
\4>k) = ET=i Uki\*i) = E*=i Ea.b UkrW ab fUt\n). From 
Lemma |^ we obtain that for each fc, r](w'(k)) = 0, where 
w'(k) a b = 52j=i UkiW l ab . The matrices Uki must therefore 
fulfill for every fc 

r r 

£ £ e ^ <bW 3 cdUkiUk] = C tJ U kt U kj = 0. (14) 

abcdi,j — l i,j*=l 

On the other hand, from Eq. ( |l3| ) we obtain 

r 

»V,//. ; 'V (15) 

fc=i 

The Slater rank 1 is thus equivalent to the existence of 
the r' x r martix U kl that fulfills Eqs. @ and 
It is convenient to represent the rows of the matrix U k i 
as vectors \R k ) in an r dimensional Hilbert space Ti. aux . 

Eqs. (p)) and (||) then reduce to YX \Rk){Rk\ = 1, 
and (R* k \C\R k ) — for all fc. One can always change 
the basis in Tl auxi i.e. replace \R k ) —* U\R k ). Such a 
transformation does not affect Eq. (|l5|), and transforms 
C — ► U T CU . Since C is symmetric, U can be choosen in 
such a way that U T CU is diagonal, and Eq. ( |l4| ) reads 
then Ei=i Ci U ki = 0. In this new basis the construction 
of U k i using the method of Wootters |Q can be carried 
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over. One can always assume that c-JJ^ is real and pos- 
itive, by chosing the phases of \Rk)- Then, one observes 
that, provided Eq. (14) is fulfilled, 



k / 



>\ci\\ul\-J2h\\ u u\ 



(16) 



Summing the above inequality over k and using Eq. (15) 
we obtain the necessary condition 



< 



E 

i=2 



(17) 



To show that it is also a sufficient condition, we take 
r' = 2 if r = 2, r' = 4 if r = 3, 4, r' = 8 if r = 5J3, and 
Uki = ±lfci exp(i8i)/ 
then all equivalent to 



The equations in Eq. (14) are 



E 

i=2 



Ci exp(2i^ 



(18) 



and the angles 6i can indeed be choosen to assure that 
Eq. (18) is fulfilled, provided the condition ( p~7|) holds. 
The Uki signs are designed in such a way that Eq. ( |l5| ) is 
fulfilled. Thus for r> = 2 we take (++), (+-) for i = 1, 2, 

for r' = 4 we take (++++), (-H ), (+-+-), (H h) 

for i = 1, . . . , 4 (or any 3 of them for i = 1, . . . , 3), and 
finally for r> = 8, (+ + + + + + ++), (+ + + H 

),(+ + -- + + ),(+ + ++),(+- + - + - 

H — ), (H 1 1 h). In the latter case we take again 

as many vectors as we need, i.e. i = 1, . . . , 5 < r < 6. 

The above Theorem is an analog of the Peres- 
Horodecki-Wootters result for two-fermion systems hav- 
ing a single-particle space of dimension 2K < 4. The sit- 
uation is much more complicated, when we go to K > 2; 
this is similar to the case of the separability problem in 
C M ® C N with MN > 6. These issues are investigated 
in Section |v|. In the following section, however, we shall 
concentrate on the case K = 2. 



we will denote such unnormalized states just as states 
occurring in the previous sections, i.e. without the bar. 

The measure (|1]) has all desired properties Jl9| , p3]| , 
such that it vanishes iff \tp) has Slater rank 1, and it is 
invariant with respect to local bilateral unitary opera- 
tions, or, in another words, with respect to changes of 
the basis in the single particle space. 

Having defined the measure for the pure states, we can 
consider the following definition: 

Definition 2 Consider a density matrix p acting in 
A(C A ® C 4 ), and all its possible convex decomposi- 
tions in terms of pure states, i.e. p = £\ \ipi)(ipi\ = 
^2iPi\tpi}{tpi\, where the unnormalized states \ipi) = 
s/Pi\t^i}] the Slater correlation measure of p, Sl(p) is de- 
fined as 



Sl(p)=inf\J2PiVm) 



where the infimum is taken over all decompositions. 

In other words, Sl(p) is the minimal amount of Slater 
correlations of the pure states that are needed in order 
to construct p, and there is a construction of p that uses 
pure states with "averaged" Slater correlation Sl(p). 

Note that SiK^dVO) = Y^iVd^i))- As we shall see 
below, the measure Sl(p) can be related directly to the 
matrix CV, in Eq. ([ll]), and to its "concurrence". It is 
invariant not only with respect to local bilateral unitary 
operations, but it also cannot increase under local bilat- 
eral operations. These are trace preserving maps of the 

where each A« 
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form p -> M (p) = J2j A j ® A jP A 

acts in C 4 , and A \ A i> ® A j A j = ^- Such transforma- 
tions correspond to mixtures of density matrices obtained 
after nonunitary changes of the basis in the single-particle 
space. It is easy to see that 



Sl(M(p)) = J2 \&*Aj\S«f) < SI (P) 



IV. SLATER CORRELATION MEASURE 

The similarity of our approach to that of Wootters jlS| ] 
can be pushed further, and in particular allows us to de- 
fine and calculate, for the case of K — 2, the "Slater 
formation measure" (in analogy to entanglement forma- 
tion measure fill). 

To this aim we first consider a pure (normalized) state 
\ip) = J2 a f, w abfafb\ty » an< ^ define the Slater correlation 
measure of \tb) as in Lemma | (cf. Ref. Q), 



*7(|Vi» = |<V#>l- 



(19) 



with being the dual of Obviously, the notion of 
dual states, as well as the function rj(.) in Eq. (|l9|), can 
be defined also for unnormalized states. In the following 



We have the following theorem: 
Theorem 2 For any p acting in A(C A (8>C 4 ) 



^(p) = M-EN> 



(20) 



where Ci are the diagonal elements of C (Eq. ) in the 
basis that diagonalizes it. 

Proof: The proof is essentially the same as the one in 
the previous section. Let us consider an arbitrary expan- 
sion of a given density matrix, p = Y^k=i \4>k)(<l>k\, where 
Ifik) — Y^=i UkjY^j)- Here denote the usual "sub- 
normalized" eigenvectors of p with | ) being equal 
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to the j-th nonzero eigenvalue of p |18|. It is easy to see 
that 



CijUkiUk-, 



(21) 



and ELi u* kl u kJ 



By changing the basis to the 



one in which C is diagonal we get (after choosing the 
phases of Uki such that c\U kl are real and positive): 



fc=i 



fe=i 



>ki-Ei 



(22) 



This inequality becomes an equality when we use the 
same construction of Ukj as in previous section, namely 
Ukj = ilfej exp(i9j)/ Vr*, with 9j selected in such a way 
that (independently of fc), 



E 



cM 3 



i 



K'i 



E 

i=2 



(23) 



The above construction provides, to our knowledge, a 
rare example of an analog of the entanglement forma- 
tion measure that can be evaluated analytically. Ob- 
viously, since we have introduced the concept of Slater 
coefficients, we may define other Slater correlations mea- 
sures for pure states in terms of appropriately designed 
convex functions of the Slater coefficients (in analogy to 
entanglement monotones [Q). For K = 2, and most 
probably only for K = 2, all those measures are equiva- 
lent and the corresponding induced measures for mixed 
states can be calculated analytically. 



V. SLATER WITNESSES 

We now investigate fermion systems with single- 
particle Hilbert spaces of dimension 2K > 4. In this case, 
a full and explicit characterization of pure and mixed 
state quantum correlations, such as given above for the 
two- fermion system with K = 2, is apparently not pos- 
sible. Therefore one has to formulate other methods to 
investigate the Slater number of a given state. We can, 
however, follow here the lines of the papers that we have 
written on entanglement witnesses JTo| , ^6[ , and Schmidt 
number witnesses |l7| ]. 

In order to determine the Slater number of a density 
matrix p we note that due to the fact that the sets Slk are 
convex and compact, any density matrix of class fc can 
be decomposed into a convex combination of a density 
matrix of class fc — 1 , and a remainder 8 p5[ : 

Proposition 1 Any state of class fc, pk, can be written 
as a convex combination of a density matrix of class fc — 1 
and a so-called fc— edge state 8: 



P k = (l-p)pk-i+p6, 1>P>0, (24) 

where the edge state 8 has Slater number > fc. 

The decomposition ( p4]) is obtained by subtracting pro- 
jectors onto pure states of Slater rank smaller than fc, 
P = \jj) <k ){i) <k \ such that p k - XP > 0. Here |-0 <fe ) 
stands for pure states of Slater rank r < fc. Denoting 
by K(p), R(p), and r(p) the kernel, range, and rank of 
p, respectively, we observe that p' oc p — \\4> <k }{ip <k \ is 

non negative iff \ip <k ) £ R(p) and A < (?/> < ' c |p~ 1 |V' < ' c ) 1 
(see ^5[). The idea behind this decomposition is that the 
edge state 8 which has generically lower rank contains all 
the information concerning the Slater number fc of the 
density matrix pk- 

As in the case of Schmidt number, there is an opti- 
mal decomposition of the form ( pi]) with p minimal. Al- 
ternatively, restricting ourselves to decompositions pk — 
J^iPil^i')^^ I with all ri < fc, we can always find a de- 
composition of the form ( ]24| ) with 8 € Slk- We define 
below more precisely what an edge state is. 

Definition 3 A fc-edge state 8 is a state such that 8 — 
e\tp <k )(ip <k \ is not positive, for any e > and \ip )■ 

Criterion 1 A mixed state 8 is a fc-edge state iff there 
exists no \ip <k ) such that \i) <k ) G R{8). 

Now we are in the position of defining a fc-Slater wit- 
ness (fc-SlW, fc > 2): 

Definition 4 A hcrmitian operator W is a Slater witness 
(S1W) of class fc iff Tr(WV) > for all a E Sl k -u and 
there exists at least one p £ Slk such that Ti(Wp) < 0. 

It is straightforward to see that every S1W that detects 
p given by ( pi| ) also detects the edge state S, since if 
Tr(Wp) < then necessarily Ti(WS) < 0, too. Thus, the 
knowledge of all SlW's of fc-edge states fully characterizes 
all p e Slk ■ Below, we show how to construct for any edge 
state a S1W which detects it. Most of the technical proofs 
used to construct and optimize Slater witnesses are very 
similar to those presented in Ref. (lO) for entanglement 
witnesses. 

All the operators we consider below act in A(C 2K (8 
C 2K ). Let 8 be a fc-edge state, C an arbitrary posi- 
tive operator such that Tr(<5C) > 0, and P a positive 
operator whose range fulfills R(P) — K(8). We de- 
fine e = inf|_^ <fc ^ (7/i <fe |P|'j/; < ' i -') and c = sup (ip\C\ip). 

Note that c > by construction and e > 0, because 
R(P) = K(8) and therefore, since R{8) does not contain 
any \ip <k ) by the definition of edge state, K(P) cannot 
contain any \4> <k ) either. This implies: 



Lemma 4 Given a fc-edge state 6, then 



W = P--C 

c 



(25) 



is a fc-SlW which detects 8. 



G 



The simplest choice of P and C consists in taking projec- 
tions onto K{5) and the identity operator on the asym- 
metric space l a , respectively. As we will see below, this 
choice provides us with a canonical form of a fc-SIW. 

Proposition 2 Any Slater witness can be written it the 
canonical form: 



W = W - el„ 



(26) 



such that R(W) — K(5), where S is a fc-edge state and 
0<e<M meSk _MWm. 

Proof: Assume W is an arbitrary fc-SIW such that 
Tr(Wa) > for all a £ Slk-i, and there 3 at least one 
p such that Tr (Wp) < 0. W has at least one negative 
eigenvalue. Construct W + el a = W, such that W is a 
positive operator on A(C 2K ® C 2K ) but does not have a 
full rank, K(W) ^ (by continuity this construction is 
always possible). But (ip <k \W\tp <k ) > e > since W is 
a fc-SIW, ergo no \^ <k ) € K{W).U 



Definition 5 A fc-Slater witness W is tangent to Slk-i 
at p if 3 a state p € <S7fc_i such that Tr(Wp) = 0. 



Observation 1 The state p is of Slater class k — 1 iff for 
all fc-SW's tangent to Sl k -i, Tr (Wp) > 0. 

Proof (See [|o|): (only if) Suppose that p is of class 
fc. From the Hahn-Banach theorem it follows that there 
exists a fc-SIW, W, that detects it. We can subtract el a 
from W, making W — el a tangent to Slk-i at some a, 
but then Tv(p(W - el)) < 0.D 



can optimize the witness by subtracting from it a posi- 
tive operator P, such that PT W = 0. For example, for 
Slater witnesses of class 2 this is possible provided that 
infi e \ eC 2K [Pe W e Pe ]mm > 0. Here for any X act- 
ing on A{C 2K <g> C 2K ) we define 



X K = 



[e,.\X\e, 
,e\X\e,: 



(e,.\X\.,e 
„e\X\.,e) 



(28) 



as an operator acting in C 2K , and [AT] m i n denotes its min- 
imal eigenvalue (see fTpf ). An example of an optimal 
witness of Slater number k in A(C 2K ® C 2K ) is given by 



W 



K 



(29) 



where V is a projector onto a "maximally correlated 
state", |*) = ^^/t^/t^fi) (cf . Eq .|) The 
reader can easily check that the above witness opera- 
tor has mean value zero in the states f r-\f r-\\Q) for 
i = 1,2, but also for all states of the form where 



„t _ ft Jtpxt , ft 

#1 — ■'a 1 (l) e + •'ai(2) 

i ft 1^21 i ft J' 

„t _ _ f t p-t^ia — I — ft 

f2 — J ai (l) e + •'a 1 (2) ( 



Ja 2 (l) e +/a 2 (2) e 



(30) 



(31) 



for arbitrary (fij, i,j — 1,2. The set Tyy spans in this 
case the whole Hilbert space A(C 2K ® C 2K ), ergo W is 
optimal. 



A. Optimal Slater witnesses 

We will now discuss the optimization of a Slater wit- 
ness. As proposed in p0| and |l7j] an entanglement wit- 
ness (Schmidt witness)W is optimal if there exists no 
other witness that detects more states than it. The same 
definition can be applied to Slater witnesses. We say that 
a k~ Slater witness Wi is finer than a k— Slater witness 
Wi , if Wi detects more states than Wi . Analogously, we 
define a k— Slater witness W to be optimal when there 
exists no finer witness than itself. Let us define the set 
of \ip <k ) pure states of Slater rank k — 1 for which the 
expectation value of the /c-Slater witness W vanishes: 

T w = W <fe ) s. t. {iP <k \W\^P <k ) = 0} , (27) 

i.e. the set of pure tangent states of Slater rank < k. 
W is an optimal fc-SIW iff W - eP is not a fc-SIW, 
for any positive operator P. If the set Tw spans the 
whole Hilbert space A(C 2K <g> C 2K ), then W is an opti- 
mal fc-SIW. If T w does not span A(C 2K <g> C 2K ), then we 



B. Slater witnesses and positive maps 

It is interesting to consider linear maps associated with 
Slater witnesses via the Jamiolkowski isomorphism pi[ . 
Such maps employ W acting in Ha ® Ti-s = C 2K ® C 27 ^, 
and transform a state p acting in Ha <8 He = C 2K C 2K 
into another state acting in Hb ® He = C 2K <8> C 2K , 
M(p) — Tta(WPa )• Obviously, such maps are positive 
on separable states: When p is separable, then for any 
|*) S Hb®H c , the mean value of (^\M(p)\^), becomes 
a convex sum of mean values of W in some product states 
|e, /) G Ha 8> Hb- Since W acts in fact in the anti- 
symmetric space, we can antisymmetrize these states, i.e. 
K /) — * (K /) ~ l/j e ))- Such antisymmetric states have, 
however, Slater rank 1, and all S1W of class k > 2 have 
thus positive mean value in those states. This class of 
positive maps is quite different from the ones considered 
in Refs. they provide thus an interesting class of 

necessary separability conditions. The map associated 
with the witness ( |29| ) is, however, decomposable, i.e. it 
is a sum of a completely positive map and another com- 
pletely positive map composed with transposition. This 
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follows from the fact that the witness operator has a posi- 
tive partial transpose, i.e. it can be presented as a partial 
transpose of a positive operator. 



knowledges partial support from the Swiss National Sci- 
ence Foundation. 



VI. CONCLUSIONS AND OUTLOOK 

Summarizing, we have presented a general characteri- 
zation of quantum correlated states in two-fermion sys- 
tems with a 2A"-dimensional single-particle space. This 
goal has been achieved by introducing the concepts of 
Slater deomposition and rank for pure states, and Slater 
number for mixed states. In particular, for the important 
case K — 2 the quantum correlations in mixed states can 
be characterized completely in analogy to Wootters' re- 
sult for separated qubits Jig] , and using the findings of 
Ref. (ll| for pure states. Similarly to the case of sepa- 
rated systems, the situation for K > 2 is more compli- 
cated. Therefore, we have also introduced witnesses of 
Slater number k 7 and presented the methods of optimiz- 
ing them. 

Possible directions for future work include generaliza- 
tions of the present results to more than two fermions, 
and the development of an analogous theory for indis- 
tinguishable bosons. For this purpose a lot of the con- 
cepts developed so far are expected to be useful there 
as well. However, there are certainly also fundamental 
differences between quantum correlations in bosonic and 
fermionic systems. As an example consider the notion of 
unextendible product bases introduced recently for sep- 
arated systems [M. These are sets of product states 
spanning a subspace of the Hilbert space whose orthogo- 
nal complement does not contain any product states. All 
such unextendible product bases constructed so far in- 
volve product states of the form \tfj)(S> |%) with and |x) 
being non-orthogonal. In the analogous fermionic state 
non-orthogonal contributions are obviously cancelled out 
by antisymmetrization, unlike the bosonic case. In fact, 
all explicit constructions of unextendible product bases 
known so far [ p6[ can be taken over directly to bosonic 
systems to give "unextendible Slater permanent bases" . 
These are sets of symmetrized product states spanning 
a subspace of the symmetrized Hilbert space, whose or- 
thogonal complement does not contain any such states. 
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APPENDIX A: 

We now list further properties of the correlation mea- 
sure r\ for pure states |*) = J3 a ,6=i w abfXfi of two 
fermions in a four-dimensional single-particle space [^3| , 
and add some further remarks. 

The matrix w transforms under a unitary transforma- 
tion of the one-particle space, 



as 



w i ► Uwll , 



(Al) 



(A2) 



where U T is the transpose (not the adjoint) of U. Un- 
der such a transformation, |*) i— > |$) = U\^f), scalar 
products of the form (^i]^) remain unchanged up to a 
phase, 



(#i|# a > =det I # 2 > 



(A3) 



Therefore, in particular, ry(|^)) is invariant under arbi- 
trary single-particle transformations. 

The dualisation of a state |^) can be identified as a 
particle- hole-transformation, 

U p - h f+U+_ h = f a , U p ^ h \fl) = ft ft ft ft P), 

(A4) 

along with a complex conjugation. In fact, the operator 
of dualization T>, |<3>) i— > IvP) = T>\^), can be written as 



V = —Up-hIC . 



(A5) 



where K. is the usual operator of complex conjugation 
which acts on a general state vector as 



K, (a\a) + b\/3)) = a*K\a) + b*K,[ 



(A6) 



Its action on the single-particle basis states and the 
fermionic vacuum is given by 



£ft£ = ft 



£falC = fa , £|n) = |fi). (A7) 



The relations ( |A7| ) are to be seen as a part of the defini- 
tion of K and refer explicitly to a certain single-particle 
basis defined by the operators f a , ft- However, switch- 
ing to a different com plex conjugation operator /C', fulfill- 
ing the relations ( A7) in a different basis, has only trivial 
effects without any physical significance. In particular, 
as one can see from the properties given above, the cor- 
relation measure 7y( | \J7) ) = |(^|^)|, = remains 
invariant under such an operation. 
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Eq. (A3) implies that T> is unchanged by unitary single- 
particle operations, 



UVU + =V & [U,V]=0 
which can also be expressed as 



(A8) 



(A9) 



for any unitary single-particle transformation U. 

The dualisation operator T> is the antiunitary imple- 
mentation of the particle-hole-transformation. We note 
that the complex conjugation involved there is necessary 
for T> being compatible with single-particle transforma- 
tions U, 



vuf+u+v- 1 = J2 UUb 



(A10) 



If the complex conjugation would be left out, hi and T> 
would not commute. 

The relation of the correlation measure r\ to an antiu- 
nitary operator is similar to Wootters' construction for 
a separate system of two qubits Q. The correlation 
measure there ("concurrence") relies on the time inver- 
sion operation. The operator of time inversion in the 
two-qubit system is invariant under local unitary trans- 
formations in each qubit space. This property is similar 
to the invariance of the dualisation operator under uni- 
tary transformations in the single-particle space. 
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